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ON THE ESTIMATION OF RELATIVE HUMIDITY PROFILES 
FROM MEDIUM RESOLUTION INFRARED SPECTRA 
OBTAINED FROM A SATELLITE 

Barney J. Conrath 
ABSTRACT 

Measurements of radiation emitted by the atmosphere of the earth in an in- 
frared absorption band of water vapor, obtained with satellite borne instrumenta- 
tion, contain information on atmospheric relative humidity. Two methods are 
developed for estimating tropospheric relative humidity profiles from infrared 
spectral measurements for which the spectral resolution elements are narrow 
compared to the total width of the absorption band, but wide compared to a single 
absorption line. The methods, which are essentially compelmentary, consist of 
a direct estimation technique which requires a minimum of a priori knowledge of 
the behavior of the relative humidity profile and a statistical estimation technique 
which can make full use of a knowledge of the statistics of tropospheric humidi- 
ties in situations where such knowledge exists. An analysis of the propagation 
of errors in the measured spectral intensities indicates that meaningful esti- 
mates should be obtained from the 6.3 micron water vapor band in the presence 
of realistic instrumental noise for most types of atmospheres. One exception is 
the polar winter atmosphere where catastrophic error propagation occurs be- 
cause of the behavior of the temperature profile. An examination of the effects 
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on the inferred relative humidity profile of errors in the temperature profile 
employed in the estimation reveals that the relative humidities inferred in the 
lowest layers of the troposphere are extremely sensitive to errors in the surface 
temperature, and this may prove to be the limiting factor in obtaining complete 
relative humidity profiles. Examples of applications of the techniques to syn- 
thetic data from model atmospheres and to real data obtained with a balloon 
borne infrared interferometer spectrometer are given. 
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ON THE ESTIMATION OF RELATIVE HUMIDITY PROFILES 
FROM MEDIUM RESOLUTION INFRARED SPECTRA 
OBTAINED FROM A SATELLITE 


I. Introduction 

Remote radiometric measurements of the thermally emitted infrared radia- 
tion from the atmosphere and the surface of the earth, performed with satellite- 
borne instruments, contain information on a number of parameters of interest in 
atmospheric physics and meteorology. Among these parameters are the tem- 
perature profile and the vertical distribution of non-uniformly mixed optically 
active gases. A considerable literature exists on the problem of estimating 
temperature profiles from radiometric measurements (see Wark and Fleming, 
1966, and references therein), and the problem of obtaining the vertical distri- 
bution of non-uniformly mixed gases was treated formally by King (1963). In the 
present study, we shall be concerned with the specific problem of obtaining esti- 
mates of the vertical humidity distribution in the troposphere. 

Several in the TIROS series of meteorological satellites carried instrumen- 
tation for making measurements in the 6 to 6.5 micron water vapor absorption 
region and the 8 to 12 micron "window" in addition to other passbands. These 
measurements have been used in estimating average tropospheric relative 
humidities, employing a method developed by Moller ( Moller , 1961, 1962; Moller 
and Raschke, 1964; Raschke and Bandeen, 1967). Recently Smith (1967) has 
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developed a technique for making estimates of the water vapor mixing ratio as 
well as temperatures in the troposphere for a proposed 5-channel radiometer 
experiment. 

Several forthcoming meteorological satellites will carry infrared inter- 
ferometer spectrometer (IRIS) experiments (cf. Hanel and Chaney , 1965). The 
first of these will cover the spectrum from 5 to 20 microns with a resolution of 
5 cm" 1 . A spectrum obtained with an instrument of this type during a high alti- 
tude balloon flight near Palestine, Texas, on 8 May 1966 is shown in Figure 1. 
Included in this spectral region is the 6.3 micron water vapor absorption band. 
The possibility of having such measurements on a global basis provides motiva- 
tion for the development of the theory of the estimation of tropospheric relative 
humidity profiles from spectra for which the individual spectral resolution ele- 
ments are small compared to the total width of the absorption band but large 
compared to an individual absorption line. It is the purpose of the present study 
to investigate certain aspects of that theory. 

In the development of any complete inversion technique a large number of 
factors must be taken into consideration. Here we shall confine ourselves to the 
development of basic approaches to the problem. The general principles of con- 
stituent inversion will be considered, and a direct estimation method will then 
be developed and used to study the sensitivity of the relative humidity estimates 
to errors in both the intensity measurements and in the temperature profile used 
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in the inversion. A statistical estimation method will then be examined. Ex- 
amples of applications to synthetic data calculated from model atmospheres and 
real data from an IRIS balloon flight will be given. The effects produced by the 
possible presence of particulate matter will not be included in the scope of this 
study. 

II. Basic Principles 

For a nonscattering atmosphere in local thermodynamic equilibrium the 
spectral intensity at the top of the atmosphere can be written from the solution 
to the radiative transfer equation in the form 

f X s 

I O') =B(zaT s )tO,X s ) - J BkTtX)] (1) 

Xco 

where B (z/,T) is the Planck function at wavenumber v and temperature T, X is 
an arbitrary independent variable taken as increasing downward toward the 
surface, and r ( v, X) is the transmissivity at wavenumber v of the column of gas 
between level X and the top of the atmosphere. The subscript s refers to sur- 
face values. It is assumed in (1) that the surface has unit emissivity, so the 
contribution from the surface is given by the first term. The second term in 
(1) represents the atmospheric contribution and for a given wavenumber can be 
regarded as an average of the source function over the weighting function 

( v,X)/dX which is everywhere positive since r decreases with increasing X. 



% 

$ 

7 

The principle of obtaining information on the atmospheric water vapor 
content from measurements of I (v) can be seen in the following way. Since 
t(i^,X) at any level X depends on the amount of water vapor above that level, a 
change in the water vapor content of the atmosphere will cause a change in the 
weighting function relative to the temperature profile so a different segment of 
the source function will be sampled, and I (^) will be changed. This effect is 
illustrated in Figure 2 where weighting functions in the 6.3 micron water vapor 
absorption band are shown for two model atmospheres with the same temperature 
profile but different relative humidities. To obtain information at several dif- 
ferent atmospheric levels, measurements must be made at several points in the 
absorption band ranging from the more opaque band center to the less opaque 
band wings. Weighting functions at two different points in the 6.3 micron water 
vapor band are included in Figure 2. The calculations employed in Figure 2 
were made using the water vapor transmissivities given by Moller and Raschke 
(1964) and correspond to resolution elements about 40 cm -1 wide. The sampling 
of the source function by the weighting function is not very sensitive to structure 
in the humidity profile which is of a scale small compared to the characteristic 
width of the weighting functions themselves. Any attempt at retrieving this fine 
structure can result in an instability against errors in the measured intensities 
just as occurs in the temperature inversion problem. Thus, the information 
obtainable on the humidity profile will be limited by the noise in the measure- 
ments of I (v ). 
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TTT . Direct Estimation 

Throughout the remainder of this paper we shall utilize the transmissivities 
of Moller and Raschke (1964) in all of our examples. The use of these functions 
limits us to consideration of spectral resolution elements no narrower than 
about 40 cm *. However, the techniques considered should be equally applicable 
to data of higher spectral resolution when used with appropriate transmissivities. 

There are a number of approaches which can be taken to the water vapor 
inversion problem; we shall consider two of these here. The first of these 
might be called a direct estimation of the relative humidity profile. 

Let us assume that we have available to us measurements of outgoing spec- 
tral intensities, and call this measured spectrum I (^). Let us further assume 
we have available to us the temperature profile T (X). Now if we arbitrarily as- 
sume some representation for the relative humidity profile containing a number 
of free parameters , a 2 , . . . a n and substitute this into (1) we obtain a para- 
metric representation of the spectrum in terms! of the a's. For example, the a's 
may be constants in some analytic form or coefficients in an expansion in terms 
of some function set. For convenience of notation let us denote the measured and 
eoretical intensities at the ith wave number by I i and I i respectively and 
efine the n-dimensional vector a whose components are a x , a 2 , ... a n . We can 
*y to evaluate a and hence get an estimate of the relative humidity profile by 
,‘tting the measured intensities equal to the theoretical representation at m 
points (m > n ) in the absorption band; i.e., 
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I i («) 


i = 1, 2, . . . m 


( 2 ) 


Since in the practical situation we can extract only a few independent pieces of 
information with many of the measured intensities being redundant to within the 
experimental error, we will generally have many more measurements than 
parameters to be determined in experiments such as the IRIS. Hence ( 2 ) will be 
an overdetermined generally nonlinear system of algebraic equations. 


A number of techniques exist for evaluating such a system of equations. 
The method which we have chosen is a generalization of the Newton-Raphson 
method. A Taylor expansion is made about some first guess, say a 0 , and 
truncated after the linear term. This gives 


where 



2 , 


m 


AI . 

1 


r 


- I s (a 0 ) 


( 3 ) 


and 


Aa . - a. - a 0 

j j j 


An ordinary least-squares solution can now be obtained for Aa from the over- 
determined linear system (3) which can be expressed in the well known matrix 
form 

Aa - (A* A)' 1 A* Al ( 4 ) 
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where 

A.. = (Bl ./3a . ) o 

and A* is the matrix transpose of A. The corrections Aa can be used to obtain 
what is hopefully an improved approximation to a which can be used as a new 
first guess and the process iterated until some convergence criterion is satis- 
fied, assuming convergence is obtained. 

The above procedure is equivalent to making a least-squares fit of I ( v ; a ) 
to the measured spectrum I (v). This does not imply that the resulting estimate 
for the relative humidity profile is necessarily a best fit to the true profile in a 
least-squares sense. One would prefer to have a technique which provides a 
least-squares fit to the true relative humidity profile, but of course this cannot 
be done directly since the true profile is unknown. However, it can be done in a 
statistical sense, and this approach will be considered in the following section. 

The estimates obtained with this technique will generally depend on the 
representation assumed. There will be two types of errors involved, one of 
which is that due to the inability of the representation chosen to fit the profile 
exactly even for perfect measurements. The other type of error is due to 
errors in the inferred values of the parameters of the representation, resulting 
from errors in the measured values of the spectral intensities. In order to 
make the representation error as small as possible, a priori knowledge of the 
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humidity profiles can be used to make a reasonable choice of representation. 

One particular class of representations is that in which the relative humidity 
profile r (X) depends on the parameters in a linear fashion, i.e., 

n 

r(X)^Y^ 3 . (X) (5) 

i = 1 

where the <£’ s are arbitrary functions of X which might be members of some 
orthogonal function set or might be chosen such that r (X ) is expressed in terms 
of step functions or ramps. We shall consider examples using the form (5), but 
such a restriction is not a necessary one. 


The transmissivities which we are employing for our examples can be 
written (Moller and Raschke, 1964) 


tOaX) 


exp 


{- 


1.97l y u*(X) 
[1+6.57 ^u*(X)] 



( 6 ) 


where / i v is the generalized absorption coefficient at wave number v, and u* 
is the reduced absorber mass 


du* =(p/p 0 ) 0 - 72 



(7) 


with T 0 and p 0 referring to standard temperature and pressure. The temperature 
factor V T 0 /T has been neglected in the present study. It has been found con- 
venient to choose the independent variable x as the atmospheric pressure p for 
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the present application. Use of the form of the transmissivity given by (6) al- 
lows us to express in a rather simple form the quantities 31 j/’ra . required for 
each step of the Newton-Raphson calculation. It is convenient to start with the 
integration by parts form of (1) 


f r « 3B.(p) , , 

Ii = Bi(0) + j — — ^ i (P) dp 


where the subscript i denotes values for the ith spectral resolution element in 
the absorption band. Making use of the well known approximate relationship 
between specific humidity and relative humidity 


w (g/kg) - 622 — r 
P 


and relation (7), along with the expression for the absorber mass 


du = — wdp 
g 


we obtain from (8) 


iii- f 1 

3a. J 


BBjCp) 3t.(p) 

—Z — (P) dp 

3 P 3u* 


where 


IS 


'\o.72 e (T) 


- — <P j (P') dp' 


9 


13 


The gravitational acceleration is represented by g, and e s is the saturation 
vapor pressure which is a function of temperature only. The factor 3u* 
occurring in the integrand in (11) is to be evaluated at each iterative step using 
the value of the relative humidity profile estimated in the previous step. 

In order for such a computational scheme to be of value, the region of con- 
vergence must be sufficiently large to allow one to make a reasonable first 
guess. It appears that the question of convergence can best be examined 
empirically. To obtain some feeling for the behavior of this method of estima- 
tion, synthetic data were calculated and inversions performed for a number of 
model atmospheres. 

Application to a rather wet mid-latitude model is shown in Figure 3. A 
3 -parameter representation consisting of two ramps linear in pressure with the 
break point at the 500 mb level was employed. Synthetic data for nine 40 cm 1 
wide spectral resolution elements between 1200 cm" 1 and 1520 cm -1 were used 
with a "first guess" of a completely saturated atmosphere with 100% relative 
humidity at all levels. Convergence to the solution shown was obtained in about 
six iterations. For most model atmospheres tried, the a's were found to change 
less than 1% after five or six iterations. In practice, one could probably make 
better first guesses than the crude one used here. 

An inversion of synthetic data from a tropical model is shown in Figure 4. 
Once again a 3-parameter 2 -ramp representation with the break point at 500 mb 
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was employed. Because of the higher atmospheric opacity due to the larger 
amount of water vapor present, it was necessary to go further into the band wing 
in order to insure some sensitivity to the lower-most layers of the troposphere. 
Ten spectral resolution elements between 1160 cm -1 and 1520 cm were 
employed. 


One very important aspect of any inversion technique is the stability of the 
estimation against random errors in the measured intensities. The root-mean- 
square fluctuation in the estimated values of the parameters a due to the 
presence of a given rms error in the intensities o l can be approximated for 
small errors by 



The factors da . /BI . required in (13) are just the elements of the matrix 

(A*A ) -1 A* of (4) which are required in the inversion. If the cr are the same 

j 

for all spectral resolution elements we can define an error amplification factor 


such that 



a 

a . 


= K 

1 


I 


(14) 


(15) 


In the two ramp representation used in the examples above, the three parameters 
solved for were the relative humidities at the 100 mb, 500 mb, and 1000 mb levels. 
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The error amplification factors corresponding to these three parameters for the 
mid-latitude and tropical models used above are given in Table I. 

Table I 


Error Amplification Factors K. (erg cm* 2 sec" 1 ster -1 cm) 1 




Model Atmosphere 


rixL u-mcicT 

Mid-Latitude 

Tropical 

Polar 


r r 100 mb 

0.376 

0.334 

- 

2 -Ramp ■< 

r 500 mb 

0.231 

0.237 

- 


^ r i000 mb 

0.603 

0.730 

- 


f r l 00 - 5 75 mb 

— 

— 

0.102 

2 -Layer < 

1 

l r 575- 1000 mb 

- 

- 

1.246 


The error which is due to the representation alone will be dependent on how 
well the assumed form can be made to fit a given profile, so will vary from case 
to case. Some idea of typical errors of representation can be gotten from 
Figures 3 and 4. The IRIS experiment jorovides an example of the random error 
in the measurement of the spectral intensities which might be expected in the 
practical situation. The rms error expected in that experiment is about 0.5 
erg cm* 2 sec" 1 ster" 1 cm. In order to meet the constraints of the transmissivi- 
ties employed in the examples above, the IRIS data would have to be averaged over 
about eight 5 cm 1 wide resolution elements , so the resulting effective rms 
error would be about 0.2 erg cm" 2 sec* 1 ster" 1 cm. The effects of such an 

error on the three parameters of the inversion are shown by the error bars in 

12 








Figure 4. In this particular case the error due to noise in the measurements 
is about the same size as the representation error. 

Figure 5 shows an inversion of synthetic data for a polar winter model 
atmosphere which is probably a worst case for inferring relative humidities. 

The representation used in this case was a step function consisting of two layers 
of constant relative humidity divided at the 575 mb level. Once again the error 
bars indicate the effects of a 0.2 erg cm -2 sec -1 ster _1 cm random error in the 
measured intensities. It can be seen that the inferred relative humidity for the 
bottom layer is quite unstable against errors in the measurements. This is 
primarily due to the behavior of the temperature profile which is also shown 
in the figure. The average temperature in the bottom layer is very near the 
surface temperature so that layer looks very much like an extension of the 
blackbody surface and the outgoing intensities are highly insensitive to the 
actual value of the relative humidity in the layer. 

A quantity of considerable interest which can be calculated from the relative 
humidity profile along with the temperature profile is the total water vapor in 
an atmospheric column. The inferred total water vapor using exact intensities 
is compared with the actual value in each of the examples shown in Figures 3, 

4, and 5. The lack of complete agreement is due to the representation errors 

in the relative humidity profiles. The percentage rms error in the inferred total 

- 1 - 2 _ 1 

amount of water vapor due to rms intensity errors of 0.2 erg sec cm ster 
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cm is 12% for the mid-latitude model, 9% for the tropical model, and 41% for the 
polar winter model. 

In the examples considered, we have used representations for the relative 
humidity profile containing considerably fewer parameters than the number of 
spectral intensity measurements employed in the inversions. Attempts at using 
more complex representations containing additional parameters, lead to such 
strong propagation of instrumental errors that the results would be of little 
value. The instrumental noise of 0.2 erg sec 1 cm ~ 2 ster -1 cm for intensities 
averaged over 40 cm _1 constitutes an error of almost 10% in the most opaque 
parts of the 6.3 micron band. This is due primarily to the low values of the 
Planck intensities in this spectral region (see Figure 1). It is only through the 
use of a number of spectral intensities in a redundant sense such as in a least- 
squares calculation that we are able to obtain inversions in the presence of such 
errors. 

Up to this point we have assumed an exact knowledge of the temperature 
profile. In the practical situation, we will presumably have to rely upon temper- 
ature profiles estimated from spectral measurements in the 15 micron carbon 
dioxide absorption band, for example. Since these temperature estimates will 
contain errors, it is important to investigate the effects of errors in the temp- 
erature profile on the relative humidity estimates. 


& 
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It was pointed out by Moller (1962) that the intensity at the top of the 
atmosphere for a spectral interval in a water vapor absorption band is not de- 
pendent on the temperature profile provided the following conditions are satisfied: 
(a) the relative humidity is constant at all levels to which the outgoing intensity 
is sensitive, (b) the optical thickness of the atmosphere is large enough to pre- 
vent there being a significant contribution from the surface, (c) the temperature 
lapse' rate is constant, (d) the atmospheric transmissivity from a given level to 
the top of the atmosphere depends only on the total amount of water vapor between 
that level and the top of the atmosphere. When these conditions are met, the 
same temperature will always occur at the same optical depth regardless of the 
absolute value of the temperature at any given height. Conditions (a), (c), and 
(d) will never be rigorously satisfied in the practical situation, and (b) is not 
satisfied when we use data from the wings of the absorption band in an effort to 
get information on the lower-most layers of the troposphere. Therefore, an 
empirical investigation was made of the effects of perturbations in the tempera- 
ture profile on the estimated relative humidities. 

Temperature perturbations of -1° K were introduced successively into the 
region between 50 mb and 300 mb, the region between 350 mb and 550 mb, the 
region between 600 mb and 950 mb, and the surface temperature for a mid- 
latitude model atmosphere, and a relative humidity inversion was performed for 
each case as well as the case of the exact temperature profile. The perturbations 
resulting for the three relative humidities forming the parameters for the two 
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ramp representation are given in Table II. The effect on estimations of the total 
water vapor amount is also indicated. It can be seen that the relative humidity 
estimates are not particularly sensitive to small errors in the atmospheric 
temperature, but are quite sensitive to errors in the surface temperature. The 
reason for this extreme sensitivity is that in using the spectral region between 
1200 cm" 1 and 1520 cm 1 rather transparent portions of the band wings have 
been included in order to obtain information on the lower-most layers of the 
troposphere. Here the boundary term in (1) is large compared to the atmospheric 
contribution so the surface temperature is the parameter to which the intensities 
are most sensitive. This sensitivity can be reduced by not using the intensities 
as far out into the band wings; however, this will reduce the information gained 
on the humidity in the lower part of the troposphere. It should be noted that 
while reasonable estimates of relative humidity might be made when there are 
errors in the estimated atmospheric temperatures, the resulting estimates of the 
absolute humidity may contain large errors due to the strong sensitivity of the 
saturation vapor pressure to the temperature. 

If the temperature profile used in estimating the relative humidity profile 
must come from an inversion in the 15 micron carbon dioxide absorption band, 
there is an additional complication due to the finite water vapor absorption in 
this region as well as the "window" region near 11 microns from which the sur- 
face temperature must be obtained. Thus, the water vapor and temperature 
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inversions are essentially coupled. One approach to this problem is an overall 
iterative scheme which takes advantage of the relatively weak dependence of the 
window and 15 micron band intensities on the humidity profile. A first guess at 
the relative humidity profile is made and the surface temperature and atmospheric 
temperature profile are estimated from the window and 15 micron regions. These 
temperatures are then used along with measured intensities in the water vapor 
band to obtain an improved relative humidity profile. The process is then 
iterated until convergence is obtained. An example of such a calculation on syn- 
thetic data from the mid-latitude model used above is shown in Figure 6. Five 
iterations were required to obtain the solution shown. The representation as- 
sumed for the temperature profile consisted of two ramps linear in the logarithm 
of the pressure with the break points chosen at the 200 mb and the 800 mb levels 
and the temperature above 200 mb assumed constant. 

As a final example of a direct estimation of relative humidity, an inversion 
of actual IRIS data was performed. The data were obtained during a high altitude 
balloon flight 8 May 1966 near Palestine, Texas ( Chaney , et al., 1967). The re- 
sults are shown in Figure 7. 

IV. Statistical Estimation 

Up to this point we have used the a priori information available on humidi- 
ties in the troposphere only as a guide in picking reasonable representations for 
approximating the relative humidity profile. For those geographic regions for 
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which a large amount of information is available, it would seem desirable to 
make better use of this information. One possible approach would be to employ 
as the function set in expansion (5) empirical orthogonal functions ( Obukhov, 

I960; Holmstrom , 1963) constructed from radiosonde data collected in the past 
at the station in question. This method has been applied to the temperature 
inversion problem (Alishouse et al, 1967; Wark and Flemming , 1966). Another 
approach is to attempt to devise a regression relation between the measured 
intensities and the relative humidity at each level. This method will be investi- 
gated here. 

As mentioned earlier, the least squares fitting of the measured spectrum 
does not in general imply a least squares fit of the representation for the 
relative humidity profile to the true relative humidity profile. We can, how- 
ever attempt to minimize the mean square deviation of the estimated profile 
from the true profile in a statistical sense. A formulation of this type has 
been given for general linear systems by Foster (1961) and for the tem- 
perature inversion problem by Rodgers (1966) and by Strand and 
Westwater (1968). The procedure below follows the general derivation 
given by these authors. 

Let us consider an ensemble consisting of q measurements of the m- 
component intensity vector I and the corresponding n-component relative 

''V/ ^ 

humidity vector r . Define the m by q matrix A such that each column corresponds 


19 



23 


to a member of the ensemble of intensity measurements from which the ensemble 
mean has been subtracted. Similarly, let p be an n by q matrix whose columns 
correspond to the measured relative humidity profiles minus the ensemble 
mean. Now define a linear estimate p to the true relative humidity profiles of 
the form 

P - H<3 (16) 

where H is an n by m matrix of coefficients to be determined. Even though the 
relationship between relative humidities and spectral intensities is generally 
nonlinear, the assumed linear form should give reasonable results at least in the 
vicinity of the ensemble mean. If it is found necessary, higher order terms can 
be added to (16). 

The elements of H can be determined by requiring that the mean square 
deviation of the estimated profiles p from the true profiles p be a minimum 
which is equivalent to requiring that the diagonal elements of the matrix 

Q =- (p - p) G° - p)* (17) 

q 

be minimized. This results in 

H = CR' 1 (18) 

where C is the n by m cross-covariance matrix 

C -- - p\ (19) 

q 
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and R is the n by n covariance matrix of the measured intensities 

R = — 11 * (20) 

q 

To be able to calculate C , we need to know the true relative humidity profiles P 
when in fact all we have available are the measured values P which will gen- 
erally contain errors 

P = P + V (21) 

where V is the matrix of errors associated with the measurements. If the errors 
are random with zero ensemble mean, then p can simply be replaced with p in 
(19). 

Since it does not appear to be practical to assemble the necessary ensemble 
of both intensity measurements and corresponding relative humidity measure- 
ments, we must rely on the measured relative humidity profiles alone and use 
radiative transfer theory and an a priori knowledge of the instrumental noise in 
the spectral measurements to obtain C and R. Let us assume the measured 
intensities & can be represented by 

<T= <3 + e (22) 

where cl is the matrix of true intensities, and e represents the matrix of instru- 
mental noise, assumed to have zero mean and to be uncorrelated with the 
measured relative humidities. The calculation can be simplified by expanding 
the theoretical relationship between the intensities and the relative humidity 
profile about the ensemble mean and truncating after the linear term, i.e., 
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^ - Pp (23) 

where P is the m by n Jacobian matrix whose elements are of the form BI. /Br . 
and are evaluated using the ensemble mean relative humidity profile. This ap- 
proximation is consistant with (16) and allows us to write C and R explicitly in 
terms of the covariance matrix of the true relative humidity profiles. In cal- 
culating the true values of A required in (22), we need the true values of p , but 
again we have available only the ensemble of measurements p. By utilizing (21) 
and the previous assumptions on 17 , along with the additional assumption that r l 
and e are uncorrelated, we obtain 

H^SP*(PSP* + N)" 1 (24) 

where N is the covariance matrix of the instrumental noise 

N — t e * (25) 

q 

and S is the covariance matrix of the true relative humidity profiles 

s - (p P* - r) r*'> (26) 

q 

Thus, if the covariance matrix l/q 1777* is known, the effect of the random errors 
in the measured relative humidity profiles p can be taken into consideration in 
the computation of S. 

We can now use H in the estimation of a relative humidity profile from 
actual intensity measurements I , providing the ensemble used in calculating S 
is representative of the conditions to which the intensity measurements pertain. 
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The estimated relative humidity profile is then given by 

i -<r> + SP*(PSP* +N) _1 [I - I (<r))] (27) 

where (r) denotes the ensemble mean relative humidity profile, which is used 
along with the temperature profile in the computation of P and I( (r) ). The 
residual variances Q i i are given by the diagonal elements of the residual co- 
variance matrix (17) evaluated using (16) and (24), 

Q- S - SP* (PSP* +N)' 1 PS (28) 

As the previous authors have pointed out, the inversion is stable against instru- 
mental noise in the sense that o -S. . as N - 00 . 

X 11 11 

In order to obtain some feeling for the behavior of such an estimation scheme, 
we have employed an ensemble of twenty-five humidity profiles based on radio- 
sonde data taken at Guam Island during the months of July and August over a 
three year period. The covariance matrix for the relative humidities was cal- 
culated, and the ensemble mean profile along with the standard deviation for each 
level ( = /S u ) is shown in Figure 8. In many situations, the noise covariance 
matrix will reduce to a single noise variance times the unit matrix. The 
residual standard deviation v/Q. . at each level is shown in Figure 8 for the case 
where a n = 0.2 erg cm sec ster cm, demonstrating the amount by which 
the variance is reduced in this particular case by the spectral data over that 
for the a priori statistics alone. 
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Figure 9 shows the results of an estimation using synthetic data from a 
model atmosphere which employed as a relative humidity profile a member of 
the original ensemble which lies fairly close to the ensemble mean profile. The 
fit in terms of both the relative humidity profile and the total water vapor amount 
is quite satisfactory. A similar inversion on synthetic data is shown in Figure 
10. In this case the model atmosphere is based on a member of the ensemble 
which differs considerably from the ensemble mean. The fit is not quite as 
good as in the previous case, especially above the 500 mb level; this is partly 
due to the assumption of the linear relationship (16) which becomes increasingly 
poor the further the sounding is from the ensemble mean. 

To obtain some feeling for the dependence of the results of an estimation on 
the type of ensemble employed relative to the type of profile sought, the covari- 
ance matrix and ensemble mean for the tropical ensemble were used to perform 
an inversion on the synthetic data from the mid-latitude model considered pre- 
viously. The results are shown in Figure 11. The fit is not as good as in the 
previous examples, presumably because the sounding is atypical with reference 
to the ensemble used. 

V. Summary and Discussion 

Two possible approaches to the problem of estimating tropospheric relative 
humidity profiles from medium resolution infrared spectra obtainable from an 
earth satellite have been developed. The direct estimation technique should be 
applicable in those situations for which little a priori knowledge of the behavior 
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of the humidity profile is available. For applications to regions for which a 
considerable amount of information exists on the statistical behavior of the rela- 
tive humidity profile, it is highly desirable to take this information into consider- 
ation in performing inversions. The statistical estimation technique which has 
been described here provides one means of accomplishing this. Hence, the two 
methods are essentially complementary. 

Calculations using both model atmospheres and actual data from a balloon 
borne infrared interferometer spectrometer experiment indicate reasonably 
stable inversions can be obtained using the direct estimation technique with a 
three parameter representation, in the presence of realistic instrumental noise. 
One exception is the polar winter type atmosphere for which the measured in- 
tensities are insensitive to the humidities because of the behavior of the tem- 
perature profile. 

An empirical study of the influence of errors in the temperature profile on 
the relative humidity estimates indicates that the relative humidity inferred is 
not extremely sensitive to errors in the atmospheric temperatures. However, 
absolute quantities such as mixing ratios obtained from these relative humidities 
and temperatures will be quite sensitive to temperature errors because of the 
strong dependence of the saturation vapor pressure on temperature. When the 
more transparent wings of the absorption band are employed in an effort to ob- 
tain information on lower tropospheric layers, it is found that the inferred 
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relative humidities are very sensitive to errors in the surface temperature. 

This may prove to be the limiting factor in the accuracy of estimates of complete 
relative humidity profiles. The surface temperature must be determined in 
satellite experiments from measurements in an atmospheric "window" such as 
that near 11 microns. The temperatures obtained in this way will be influenced 
by instrumental errors and by uncertainties in the atmospheric transmissivities 
in the window region. The window transmissivities will be dependent not only 
on water vapor, but also on the particulate matter present which will generally 
be an unknown parameter. 

The linear statistical estimation technique when applied to synthetic data 
from model atmospheres shows sufficient promise to warrant further develop- 
ment in the future. The success of the method depends strongly on having a good 
estimate for the relative humidity covariances for the location in question and 
also the covariances for the noise associated with the instrument employed. 
Further study is required to determine those geographical regions for which the 
method might be reliably used. Other areas of future investigation include the 
possibility of obtaining an improvement by including nonlinear terms in (16) and 
the possibility of combining measurements in the 6.3 micron band, the window 
region, and the less opaque parts of the 15 micron band to perform a combined 
statistical estimation of tropospheric temperature and water vapor. 
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This study has concentrated primarily on the 6.3 micron water vapor band 
because the first of the satellite experiments to which the work is applicable 
will include that band. However, the basic techniques developed here should be 
equally applicable to the rotation water vapor band beyond 20 microns and that 
spectral region should not be neglected in future work. The transmissivities of 
Moller and Raschke have been used here, but other transmissivities such as 
those given by Williamson and Houghton (1965) could be fit into the computational 
techniques equally well. 

There are a number of factors which must be taken into consideration in 
the development of methods for obtaining relative humidity profiles from satellite 
measurements on a truly global basis which have not been included within the 
scope of this study. Perhaps foremost among these is the problem of treating 
the case in which the field-of-view of the instrument is partly cloud filled. 

Smith (1967) has proposed a technique for treating this problem based on the 
availability of additional information in the form of spatial scans within the 
principal area for which the inversion is to be performed. For experiments 
for which this additional information is not available, other techniques must be 
developed. Also, quantitative estimates of the effects of atmospheric turbidity 
which may be enhanced in spectral regions of strong water vapor absorption 
(Deirmendjian, 1960) are needed. However, the computational methods con- 
sidered here should provide a basis for utilizing the forthcoming satellite ex- 
periments to develop techniques which can eventually be incorporated into 


operational systems. 
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Figure 1— Spectrum obtained with an infrared interferometer spectrometer (IRIS) flown in a high 
altitude balloon 8 May 1966 near Palestine, Texas. 
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Figure 3-Inversion of synthetic data from a mid-latitude model 
atmosphere. A 3-parameter representation consisting of two 
ramps linear in pressure was employed. 
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Figure 4-A 2-ramp inversion of synthetic data from a tropical model 
atmosphere. The error bars indicate the rms dispersion in the in- 
ferred profile due to noise in the intensities with an rms value ty pi 
cal of that expected from satellite borne instruments. 
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Figure 5-A 2-layer inversion of synthetic data from a model polar 
winter atmosphere. The error bars have the same significance as 
in Fig- 4. The temperature profile employed in the model is also 
shown. 
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Figure 7-Inversion of data from the IRIS balloon flight shown in Fig 1. 
Data from a nearby radiosonde station are shown for comparison. 
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Figure 9-Stati stical estimation of a relative humidity profile using 
synthetic data from a model atmosphere. A member of the ensemble 
which did not depart greatly from the ensemble mean was employed 
as the sounding. 
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Figure 10 Stati stical estimation of a relative humidi ty profile using 
synthetic data from a model atmosphere. A member of the ensemble 
showing considerable deviation from the ensemble mean was em- 
ployed as the sounding. 
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Figure 1 1 -Stati stical estimation of a relative humidity profile using 
synthetic data from a model atmosphere. The mid-latitude sounding 
employed was not a member of the ensemble on which the statistical 
parameters of the estimation were based. 


41 


